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INTRODUCTION 



Unconventional superconductivity in cuprates is often 
measured by deviations from Bardeen, Cooper and Schri- 
effer's (BCS) phenomenology [l|. A case in point is the 
order parameter 



$(T) = £d(r,)(4 T 



(1) 



where d(rf) is the intra-layer pairing function with d-wave 
symmetry, and uniformity is assumed in suppressing the 
r dependence of <I>. In BCS theory, the order parame- 
ter is inextricably related to a gap in the quasiparticle 
excitations, whose maximal value is given by 




Ascs(T) = V$(T). 



(2) 



where V is an interaction parameter. Abcs is the pair 
breaking energy which sets the scale of the transition 
temperature T c . However, BCS theory is a mean field 
approximation which neglects all phase fluctuations. 
In underdoped cuprates, there is compelling evidence 
that T c is driven by phase fluctuations J2[ . Uemura's 
empirical scaling law T c oc Pg b (T = 0) Q and the ob- 
servation of a superfluid density jump in ultrathin un- 
derdoped cuprate films [3, @, 0] are consistent with the 
behavior of a bosonic superfluid, captured by an effective 
xy model. 

In this paper we calculate the temperature dependent or- 
der parameter of an effective Hamiltonian of charged lat- 
tice bosons (CLB). The CLB model incorporates essential 
ingredients of underdoped cuprates including extremely 
weak interlayer coupling, and long range Coulomb inter- 
actions. 

Our main result is that <!>(T) exhibits a trapezoidal shape 
in the weak interlayer coupling limit, as depicted in Fig. 
[TJ At low temperatures $(T) decreases slowly due to 
effects of anisotropic plasma frequency gaps. The effects 
of long range charge interactions, however, do not drive 



FIG. 1: Temperature dependences of normalized supercon- 
ducting order parameters. A trapezoidal shape is obtained for 
Charged Lattice Bosons (black color online) model (Eq|5jl for 
anisotropy ratio a = 10 -6 and k = 150. We see that coulomb 
interactions suppress thermal phase fluctuations relative to 
the classical xy model, Eq.©, (red color online), depicted for 
a — 10 -6 . The rapid fall toward T c is calculated within in- 
terlayer mean held theory (see text). BCS theory for d and 
s-wave order parameters is depicted for comparison, (green 
and blue colors online, respectively). 



the transition. The transition is driven by proliferation 
of vortex loops above the two dimensional Bcrczinskii- 
Kosterlitz-Thouless (BKT) Q temperature Tbkt, where 
the order paramater falls rapidly toward T c . 
Phase fluctuation theories have been previously applied 
to cuprates, with special attention to the intra-layer su- 
perfluid density pf(T) [U E3, EH . In order to explain the 
linearly decreasing temperature dependence, additional 
gapless (nodal) fermionic excitations were argued to be 
essential fl2| . 

<&(T), however, behaves differently than p^ b (T). In the 
two dimensional limit, for example, $ must vanish at all 
T > by Mermin and Wagner theorem [l3|, while p" b 



2 



jumps to a finite value below T c . Also, nodal fcrmions 
have a small effect on $(T). This is demonstrated in 
FigJU which shows s and d wave order parameters be- 
having very similarly within BCS theory. 
We propose experimental probes for the order parame- 
ter, without relying on BCS theory and Eq. ([2|). At 
weak interlayer coupling, we argue that $(T) should be 
proportional to the square root of the c-axis supcrfluid 
density. 

Angular Resolved Photocmmission Sprectroscopy 
(ARPES) finds a "pseudogap" A pg in the electronic 
spectrum, which persists well above T c 0, [HI, fl^ ]. 
Apparently, A pg (T) is not proportional to ^(T), (the 
latter of course vanishes at T c ), which violates Eq. ((2]). 
Pseudogap phenomena are often interpreted as short 
range pairing correlations well above T c . 
To address ARPES data, we employ a Boson-Fermion 
QBF) model which was derived from the Hubbard model 
jl7j | by contractor renormalization. The model describes 
the CLB system, Andreev-couplcd to fermion quasipar- 
ticles which occupy small hole pockets (or 'arcs'). Sim- 
ilar BF models were arrived at by other approaches 
[II, [01 111- Within our model, $(T) is proportional to 
the transverse nodal velocity v±(T). In Fig. we find 
reasonable agreement between the theoretical curves and 
BSCCO ARPES data for v±{T) of Refs. [H|,[22|. Fur- 
ther tests of the trapezoidal shape closer to T c would be 
desirable. 

The paper is organized as follows: The CLB model is 
introduced in Section [TTJ The order parameter is calcu- 
lated within the harmonic phase fluctuations approxima- 
tion to obtain the low temperature regime. In Section 
IIII1 the interlayer mean field theory is applied to com- 
pute the suppression of the order parameter near T c . In 
Section HVl we relate the model parameters to experimen- 
tal data for several commonly studied cuprates. Section 
IVl compares the theory to experiments, using an effective 
Boson-Fermion model to interpret the ARPES data. We 
conclude with a brief summary and discussion. 
In Appendix [A] we provide details of the analytical fit to 
the harmonic phase fluctuations result. In Appendix IBI 
we estimate the temperature region near T c where three- 
dimensional critical fluctuations are important (Ginzburg 
criterion) . 



in effect a coarse grained parameter chosen to be larger 
than the coherence length £. J is the bare intra-layer 
superfluid density, and a <C 1 is the anisotropy ratio. 
The order parameter is defined by 



*(T) = $ (cos(^ r )) 



(4) 



where $o is the zero temperature value. The quantum 
CLB model is given by 



ticlb = HxyM + ^ ^2 V ( T ~ r ') n r n r> (5) 
r,r' 

where n r is the occupation number of a charge 2e boson 
on site r, obeying the commutation relation, 



[n r , ip r i] = iS r y 



(6) 



Long range Coulomb interactions V(r) are given by the 
Fourier components 



V a 



V(r) 



167re 2 
ve b q 2 



(7) 



where v = a c is a unit cell volume and £h(q, cj q ) is the ef- 
fective dielectric function in the appropriate wave vector 
and frequency scale. 

At low temperatures, we can expand the CLB action to 
quadratic order and obtain the harmonic phase fluctua- 
tions (HPF) action, 



5, 



hpfl 



<4 + ^ (q) 



(8) 



where w„ = 2i:nT/h arc bosonic Matsubara frequencies. 
The plasmon dispersion, as derived by Kwon et. al. [Io| , 
is 



II. CHARGED LATTICE BOSONS 

The xy Hamiltonian is a lattice model of boson phase 
fluctuations, 



cos(<£ r - ip r+v 



COs(^ r - (^r+c) 



r.c 



(3) 

where r resides on a layered tetragonal lattice. 77 and 
c are in-plane and interlayer nearest neighbor vectors of 
lengths a and c respectively. The lattice constant a is 



.2/ x _ ^Ibllb+^lll 



<4(<0 



167re 2 J 

EbftPc 

l6ire 2 caJ 
ebh 2 a 2 



(9) 



where q a b and q c are the planar and c-axis wave- vectors 
respectively. 

The HPF order parameter is given by 



(10) 



3 



where the local phase fluctuations are given by, 

f d 3 g Vq / sinh(^ p (q)/r) \ 
" V J (2tt) 3 7kJ p (q) \cosh{hu p (q)/T)-l) 

(11) 

At extremely low temperatures, T -C fi^> c , all thermal 
phase fluctuations are frozen out. However, as we shall 
show in Section [Vl the experimentally interesting regime 
of large anisotropy, has a wide separation of plasma en- 
ergy scales, such that 



Variational dctrminition of h yields the IMFT equation 

h = 2a J(cos ip T ) = 2a J$2d(T, h) (16) 

where the magnetization of a single two dimensional 
layer, $>2d{T, h), is, in principle, the exact field dependent 
order parameter of the single layer CLB model. Solving 
Eq. (|16[) for h(T), yields the three dimensional tempera- 
ture dependent order parameter 

<f> im ft(T) = $2d(T,h(T)). (17) 
The transition temperature T c is given by 



fi^c < Tbkt ~ T c < hw ab (12) 

For our regime, we fit Eq. (fTTj) by the analytical approx- 
imation (see Appendix [A|l . 



(^ 2 >T,a = (jj (ai - a 2 | ln(a)|) e -^^J^/T (13) 

For the simplified case of a = c, the coefficients are given 
by: 



T c = min{T;$,„ l/( (T) = 0} (18) 

Solution of eq. (|16[) for samll anisotropics requires precise 
determination of <f>2d(/i,T) for very weak fields /i near 
T c . This is obtained by using the asymptotic critical 
properties of the order parameter near Tbkt-, which is 
not far from T c in the small a limit. 

A. BKT critical properties 



a-L w 0.045, a 2 = -0.013 

a 3 w 0.35 (14) 

Thus, expression (fT0|) reduces to the classical result of 
Hikami and Tsuneto (HT) [H, (shown later in Eq. 
in the limit T a^h^/ujab^c- In the experimentally rele- 
vant regime, ^hp/ decreases significantly slower than the 
classical model, as demonstrated in Fig. [TJ 



The two dimensional classical xy model undergoes a BKT 
transition [H at Tbkt ~ 0.89 J [HI, [HI]. Vortex pair 
proliferation changes the phase correlation temperature 
dependence from power law to exponential decay, 

— 7](T) rp rp 

(cos(^ r -^)) ~ r p . r/i{T) i, 1 T BKT (19) 

e ' s ^ 1 1 > 1 BKT 

where at low temperatures, 



III. INTERLAYER MEAN FIELD THEORY 

The HPF action ([8]) cannot describe the order parameter 
near T c since it does not include vortex excitations. In 
the narrow regime of Tbkt < T < T c proliferation of 
widely separated two dimensional vortex pairs dramati- 
cally reduces the order parameter. 

For anisotropics of order a ~ 10~ 4 — 10~ 6 , a straighfor- 
ward numerical calculation of Eq. ^ is encumbered by 
finite size limitations. Instead, we employ the interlaycr 
mean field theory (IMFT) [24| , described by a single layer 
hamiltonian in an effective field h: 



riimftih) 

n 2d (h) 



n 2d (h) 



2aJ 



2D 



2D 



— J COs((y9 r — ifir+r]) — 2h COs(tp r ) 
rr) r 

(15) 



T 
2^J' 



T <C Tbkt 



(20) 



Above Tbkt the correlation length diverges as 



£ 2 d oc exp (f3/Vi 



D 



X2d{t) = "yexp 



ft 



(3/Vt) 

(T — Tbkt) /Tbkt 
3/2, v = 7/4 



(21) 



where the exponents (3 and v were derived by Kostcrlitz 

S3- 

In order to match the transition region to the low tem- 
perature HPF order parameter, we need to determine the 
non-universal amplitude B x of X2d(T). B x was deter- 
mined numerically. We evaluated Q 2 d{h,t) by a Monte- 
Carlo simulation with Hamiltonian (|15[) . Good conver- 
gence was achieved with 10 9 spin tilts per (T, h) point, 
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Tbkt T 



FIG. 2: Determination of B x from Monte-Carlo data. The 
different fitting functions F(B X ,T) (solid lines, red color on- 
line), are defined in Eq. (|22[) . The curve with parameter value 
B x = 0.072 is chosen as the best fit to (T - Tbkt) /Tbkt- 
The dashed (gren color online) line is the two dimensional 
order parameter in the presence of an ordering field h. 

sampling every 10 5 tilts and averaging over the last 5000 
configurations. We define a fitting function 

F{B *'^ = L(J* M {KT)/{hB x ))) (22) 
The fitting procedure which is depicted in Fig. [SJ yields 

B x ~ 0.072. (23) 

For finite interlayer coupling, the classical xy model or- 
ders at T c (a) > Tbkt- Hikami and Tsuneto [23| evalu- 
ated the order parameter for small a <C 1, and obtained 

$ C/ (T) = $ a"/( 4 - 2 "' ~ $„e~«^ |lna| . (24) 

In Fig. [TJ $ C /(T) of Eq. ([2"4"| is plotted in comparison 
to the CLB model. The classical model decreases much 
faster since it does not contain the plasma gaps in the 
thermal phase fluctuations. 
The IMFT equation for T c is 

2aJ X 2d(T c ) = 1. (25) 

Using Eq. (ST) for X2d{T) and the value (23) for B x , the 
shift of T c is 

T c -T B ^r.{^-^T BKT , (26) 

The IMFT is consistent with the renormalization group 
analysis of Hikami and Tsuneto (23j. We note, however, 



that a large vortex core energy can increase the shift of 

T c above the value given by Eq. (SB) [HHl]. 

The critical field-exponent was derived by Kosterlitz p7| 

®2d(T B KT) cx h>l\ J = 15. (27) 
Combining this result with the IMFT equation (p"6[) yields 

HT BKT ) cx a 1 ^ 6 ^ = a 1 / 14 . (28) 

Thus, by Eqs. ([26]) and (28) , the order parameter drops 
rapidly between Tbkt and T c , with an average slope of 
d$(T)/dT~ -|ln(a)| 2 . 

B. Matching at the crossover 

In the crossover region, <&2<j is given by the harmonic 
mean of the temperature and field dependent singularities 
at Tbkt- 

* M (T,,,) = *«(T)(^ (29) 

Eq. (S5) correctly captures the singularities of the vari- 
ables (t, h) at the BKT transition, ho is chosen to match 
the order parameter smoothly at Tbkt-, 

h = 2aJ<P hpf {T BK T), (30) 

IMFT, as a mean field theory cannot properly capture 
three dimensional critical exponents of the xy model. 
Nevertheless, as shown in Appendix [Bj the critical 
regime, by Ginzburg's criterion is limited to 

T c -T <T BXT /|lna| 4 , (31) 

which is difficult to resolve experimentally, in the systems 
of interest. 



C. Fermionic excitations 

The CLB model ignores effects of fermionic particle-hole 
excitations, which are clearly observed in ARPES and 
tunneling. In undcrdoped cuprates, most of their spec- 
tral weight is associated with wave vectors around the 
antinodes, ((n, 0), (0, 7r)), with energies at the pscudogap 
scale A pg 3> T c . Contribution of these excitations to de- 
pletion of the order parameter temperature is of order 
T/A pg « 1. 

Nevertheless, one might worry that low energy (nodal) 
excitations might play an important role. This has been 
shown to be the case for the temperature dependence of 
the superfluid density pf{T) 0, [ll|, HH ) . 
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However, nodal excitations are weakly coupled to the 
order parameter. Consider, for example, the BCS gap 
equation, 

X=E^5)^WWT), (32) 

where 

£k = V(ek-^) 2 + M(k)A(T)|2, (33) 

and A is the BCS coupling constant. The pair wave 
function factor |ei(k)| 2 vanishes on the nodal lines k = 
(±k,k). This suppresses contributions from the nodal 
regions to the thermal depiction of the gap. As a re- 
sult, s-wave and d-wave order parameters have very sim- 
ilar temperature dependence as shown by Won and Maki 
[lil and depicted in Fig. [T] Although here we do not 
appeal to BCS theory, this observation depends only on 
the weak coupling between nodal fermions and the order 
parameter, imposed by the pair wave function symmetry. 



CLB a = 10" 6 




FIG. 3: Comparison of CLB order parameter to square root 
of c-axis superfluid density from Ref. [4l|] . Model parameters 
are a = 10 -6 , c/a = 0.5 and k = co a b/J = 150. Data was 
taken on BSCCO with T c = 87 K. Dashed line is d-wave BCS 
energy gap, given for comparison. 



IV. EXPERIMENTAL PARAMETERS 

The cuprates exhibit very large anisotropy between in- 
plane and interlayer Josephson couplings J c and J a b, 
which can be experimentally determined by the in- 
plane and interlayer zero temperature London penetra- 
tion depths A° h and A°, 



A ab — 



A° = 



/167re 2 

16ne 2 d 
Wc 



Jal 



-Jc 



(34) 



where d and a are effective lattice constants, e is the 
electron charge and c the speed of light. The anisotropy 
ratio for cuprates is in the range, 



a 



Jc_ 

Jab 



X".d 



1(T 



10" 



(35) 



Our phenomenological assignment of J a b and J c , ne- 
glects quantum corrections which become sizable near 
the critical doping toward the insulating phase. An 
alternative measure of J a b and J c is given by rela- 
tions © and the experimental measurements of ui a b and 
Lo r . by optical and microwave conductivities (cf. Refs. 

m Hi M, M, M, HI EE M, M, E3i). Thus, the 

anisotropy parameter, a, of Eq. (|35p can be determined. 
Table U contains typical experimental values of relevant 
quantities at zero temperature (except for Q and a which 
were determined viaEqs. J A 6 1) and l|35[l respectively). In 
YBCO, BSCCO and TBCCO, the interplane distance c 
is taken as the mean value. 



V. EXPERIMENTAL PROBES OF $(T) 

In cuprates, the BCS relation, does not hold, since 
the maximal gap A pg is weakly temperature dependent 
[3 [42|, while $(T) vanishes at T c . Here we propose 
experimental probes to measure $(T)/$(0). 



A. c-axis superfluid density 

Since the zero temperature interlayer pair tunneling is 
weak, the layered system can be treated as a one di- 
mensional array of Josephson junctions. Within a vari- 
ational approximation, the order parameter can be ex- 
tracted from the temperature dependence of the c-axis 
superfluid density, 



pS(T) = ^(0)|$(T)| i 



(36) 



Indeed, as seen in Fig. [31 agreement between theoretical 
curves $(T) and values extracted from electrodynami- 
cal data of BSCCO [4l[ are quite good, except near the 
transition. 



B. ARPES 

In d-wave BCS theory the quasiparticle spectrum is given 
by Eq. (J33J). Above T c , A BC s = 0, and the full Fermi 
surface should be detected as zero energy crossings of the 
ARPES quasiparticle peaks. However, in underdoped 
cuprates as temperature is raised above T c , only finite 
Fermi arcs appear around the nodal directions. The gap 
in the anti-nodal directions A„„ survives to much higher 
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Compound 


a [A] 


c [A] 


T c [K] 


/kj ai) [eV] 


h/jJc [meV] 


n [meV] 




A c [/jm] 


a [10 -4 ] 


References 


YBa 2 Cu 3 7 -5 


3.8 


5.8 


89 


1.5-2.5 


5.6-13.6 


36.1-72.3 


0.14-0.28 


1.26-7.17 


50-5 


[30, 36, 37,43J 


Bi2Sr2CaCu2 0s+5 


5.4 


7.7 


92 


0.94-1.84 


0.23-1.4 


5.6-19.5 


0.2 


110 


0.016 


[31, 32, 38, 44] 


La 2 _5Sr5Cu04 


3.8 


6.6 


40 


0.3-3 


3.7-11.2 


13-75.4 


0.19-0.28 


2-8.5 


30-3 


[16, 33, 34, 36, 39,45] 


TlaBaaCaCuaOg+i 


3.9 


7.4 


108 


1.5 


1.2-2.6 


17.4-26.8 


0.17-0.33 


2.5-8.4 


13-4 


[35,40,46] 



TABLE I: Typical planar lattice constants, a, mean interplane distances, c, critical temperatures, T c , planar and interplane 
plasma frequencies, cj a j,, ui c , energy scales Q(a,a,c) of Eq. (|13p . magnetic field penetration depths, A a (,, A c and anisotropy 
factors, a, at zero temperature. All quantities except for fl(a, a, c) and a were obtained experimentally, while Q and a were 
obtained via Eqs. ()A6[) and (|35p respectively. Some quantities depend on doping (e.g. A a (,, A c are diminished with doping) and 
values for each compound correspond to similar dopings. 




9 arc tt/2 



FIG. 4: Boson-Fermion model for the transverse quasiparticle 
excitations below and above T c . 8 is the azimuthal coordinate 
transverse to the nodal direction. Above T c (red color online), 
vanishing of on a finite 'arc' reflects the inner edge of the 
hole pocket. The pseudogap A pg is the hole fermions energy 
at the antinodal wavevectors (n, 0), (0, tt), which has no direct 
bearing on the superconducting properties. Below T c (blue 
color online), the Andreev coupling of hole fermions to hole- 
pair bosons yields a d-wave gap with a node at 8 = 0. The 
transverse nodal velocity v±(T) is a direct measurement of 
$(T). The break in the curve at the arc edge is consistent 
with 'two gaps' phenomenology [4S( |. 

temperatures [IHHH. In contrast to A pg , the transverse 
nodal velocity v± vanishes abruptly at T c pU . |22| . Below 
T c , v± introduces a singularity in the electronic prop- 
agator, which translates to an infinite correlation length 
in real space. 

A microscopic connection between v±(T) and <J>(T) can 
be provided by an effective Boson-Fermion hamiltonian 
with small hole pockets, described below. 



C. Boson-Fermion theory 

The Boson-Fermion (BF) model, which arises by a con- 
tractor renormalization of the square lattice Hubbard 
model [l7j], describes spin half fermion holes /k, s of charge 
e, coupled to the CLB as 



k,s 

+gJ2^"d(T ~ r')/ r , T /rM + h.c. (37) 

r.r' 

The last Andreev coupling term, describes disintegration 
of hole pairs into single spin-half hole fermions. In our 
version of the BF model, the fermion and boson densities, 
measured with respect to half filling, obey 

n h + 2n b = x, (38) 

where x is the total concentration of doped holes. The 
hole dispersion £k has minima near (±7r/2, ±7r/2), and 
therefore occupy four small pockets of area fraction n/j/2. 
Above T c , the small wave vector sides of the pockets ap- 
pear as the celebrated Fermi 'arcs' j47j.The pseudogap is 
given by the quasiparticle excitation energy at the anti- 
nodal wavevectors 

A P9 = 6(^,0) - (39) 

In the superconducting phase $(T) = (cos(y)). The 
hole fermions acquire the Dirac cone dispersion near the 
nodes: 

E U = ±y/(v F (k ]] - k F )f + (2 5 4>(T)fc ± ) 2 , (40) 

that is depicted in Fig. |U Thus, the transverse velocity 
directly measures the order parameter, 

wx(T) = 2 ff $(T). (41) 

In the underdoped regime, the transverse velocity is 
smaller than the pseudogap scale A pg d'(k). This is seen 
as a 'break' in at the Fermi arcs angles, as shown in 
Fig. 2J Such behavior has been observed in ARPES [48| 
and found consistent with a 'two gaps' phenomenology. 
In Fig. [5] we compare the CLB order parameter to the 
transverse nodal velocities measured on three samples of 
BSSCO by two groups [2l|, The agreement is rea- 
sonable, although the sharp break in the curves is not 
clearly confirmed. A comparison to the <i-wave BCS ex- 
pression shows a systematic trend of all the data being 
higher than BCS theory would predict. 
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CLB - a = 10" 4 




ARPES, BSCCO 

1 • T c = 92°K - 

2 ■ T c = m°K - 

3 * T c = 80°K - 



T/T c 



FIG. 5: Comparison of CLB to transverse nodal velocity, mea- 
sured on different samples by ARPES. BSSCO samples with 
T c noted in the figure. Experimental data, inclu ding error 
bars, are (1) from Ref.JlJ], (2) and (3) from Ref. [If. The 
zero temperature normalization is chosen by the lowest tem- 
perature data points. Theoretical curves for several values of 
a are drawn, using c/a = 0.5 and k = u a b/J — 150. Dashed 
line is d-wave BCS energy gap, given for comparison. 



VI. DISCUSSION 

This paper calculated the order parameter of cuprates 
using a bosonic model of hole pairs. The model includes 
crucial features of layered cuprates: long range Coulomb 
interactions and very small anisotropy ratio. It ignores 
effects of fcrmionic particle hole excitations which are 
argued to be small for $(T). The calculation predicts a 
trapezoidal temperature dependence in the small a limit, 
which is distinct from both BCS theory and the classical 
xy model. The theoretical curves are compared to data 
where the order parameter is extracted by additional the- 
oretical assumptions: the c-axis superfluid density (using 
a variational argument) and the transverse nodal veloc- 
ity (using a BF model of small hole pockets). We have 
selectively chosen data of BSCCO where a = 10~ 6 , and 
the 'trapezoidal' temperature dependence is most pro- 
nounced. In other cuprates, with larger values of a, and 
larger vortex core energies (2|| the shift Tc — Tbkt is 
larger, and the curve should be more rounded (less trape- 
zoidal) and similar to the BCS curve. 
Additional probes to $(T) could be devised. The criti- 
cal current of Josephson junctions with a higher 
T c material might be investigated. The transverse nodal 
velocity, which we have related to $ by the BF theory, 
determines the low energy tunneling spectra and Raman 
scattering (49] . In addition, it has been theoretically re- 
lated to the linear slope of the superfluid density dp a ^ jdT 
[28|, and to thermal conductivity. 

Further comparisons to experiments are warranted. 
Their success or failure may shed light on the applicabil- 
ity of the quantum lattice bosons description of cuprates 



both below and above T c . This would help us resolve 
some of the other mysteries of the pscudogap phase. 
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APPENDIX A: FITTING PHASE 
FLUCTUATIONS 

We define q 2 = q 2 + q 2 lb , if = (uj c /uj a b) 2 = aj 2 and 
7 = c/a. For ease of numerical integration Eq. (Jll[) may 
be simplified as follows 



xy <p n e 
d 3 q Vq 



(Al) 



(2w) 3 huj p (q) \cosh(f3huj p (q)) — 1 

77T 7T 

lhWab f dz [ dr 



2nU J ~~ J (z 2 + r 2 )e(r l ^) 
oo r 



sinh [e (y, I ) /T] 
'cosh[e (77, ~) /T] — l' 



(A2) 



where the last expression was obtained in cylindrical co- 
ordinates. The dispersion is thus parametrized by 



rj 2 (z/r) 



+ (z/r) 2 



(A3) 



At extremely low temperatures, T <C fiuj c , all thermal 
phase fluctuations are frozen out. However, due to the 
large anisotropy, and poor screening, there is a wide sep- 
aration of energy scales between the interplane plasma 
gap, Hu) c , and the planar gap, huj a b and it turns out that 



(A4) 



At low temperatures, the integral in Eq. (fTT]) may be 
parametrized as 



<<^ 2 > 



AT 



(A5) 



where the energy scale, f2(o!, 7), and the coefficient 
A[a 1 7) may be parametrized by 
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Q(a, 1 )^h^—+0.26^j ^J^T C , (A6) 

and 

A(a, 7) w Ai(7) - A 2 (7) Ln(a) 

^1(7) w 0.0297 + O.OI67 2 

A 2 (j) ~ 0.24 7 -O.II7 2 . (A7) 

The low temperature magnetization, (T), is given 

by 

$h P /(T) (x e"5^ 2 > = C*(T)a^ oxp(n/T) , (A8) 
where the coefficient C (T) is given by 

C{T) = e -^ex P (n/T)_ ( A9 ) 

Notably, 7 is of order unity and the energy scale, f2(a, 7), 
in Eq. (|A5j) is proportional to the geometric average of 
the interplane and planar plasma energies. 

APPENDIX B: GINZBURG'S CRITERION FOR 
INTERLAYER MEAN FIELD THEORY 

One would like to know, in which regime can we trust the 
IMFT near the transition temperature. Here we estimate 
the critical region using the standard Ginzburg Criterion. 



At small a, we see that the magnetization only varies 
rapidly below T c , in the narrow region of width AT C given 
by Eq. (|26|) . Within that region, $ im / t (T) drops from 
^>hpf(Tj(T), as given by the harmonic phase fluctuations 
(|A8[) . to zero at T c , with a mean field behavior, 

/iT-TlV 3 1 

awt ~ *h P f(T KT ) ( —£rf) - 0=2- ( B1 ) 

Ginzburg's criterion [50l. [5l|. estimates the temperature 
region below T c , where critical 3D fluctuations become 
important and IMFT breaks down. This is where order 
parameter fluctuations averaged over a correlation vol- 
ume of size V£ = £a b £ c exceed their average, i.e. 

<(A$) 2 >v, = 5(q ^°' T) = ^ > *Lft<?). (B2) 

Using the mean field estimation of £ c ~ c(|T— T c \/T c )~~^ 
and Eq. (|B1|) . the critical regime is given by 

\T-T C \< AT 2 JT C « AT C , (B3) 

which is much smaller than the already narrow region 
of AT C , where 2D vortex pair fluctuations suppress the 
order parameter. In summary, for layered systems with 
large anisotropy, IMFT theory holds up to temperatures 
very close to T c . 
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